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Abstract 

In this paper, we study the complexity of computing the determinant of a matrix over a 
non-comniutativc algebra. In particular, we ask the question, "over which algebras, is the 
determinant easier to compute than the permanent?" Towards resolving this question, we show 
the following hardness and easiness of noncommutative determinant computation. 

• [Hardness] Computing the determinant of an n x n matrix whose entries are themselves 

2x2 matrices over a field is as hard as computing the permanent over the field. This 
extends the recent result of Arvind and Srinivasan, who proved a similar result which 
however required the entries to be of linear dimension. 

• [Easiness] Determinant of an n x n matrix whose entries are themselves d x d upper 

triangular matrices can be computed in poly(n'') time. 

Combining the above with the decomposition theorem of finite dimensional algebras (in partic- 
ular exploiting the simple structure of 2 x 2 matrix algebras), we can extend the above hardness 
and easiness statements to more general algebras as follows. Let A be a finite dimensional 

algebra over a finite field with radical R{A). 

• [Hardness] If the quotient A/R{A) is non-commutative, then computing the determinant 
over the algebra A is as hard as computing the permanent. 

• [Easiness] If the quotient A/R{A) is commutative and furthermore, R{A) has nilpotency 
index d (i.e., the smallest d such that R{A)'^ = 0), then there exists a poly(n'')-time 
algorithm that computes determinants over the algebra A. 

In particular, for any constant dimensional algebra A over a finite field, since the nilpotency 

index of R{A) is at most a constant, we have the; following dichotomy theorem; if A/R{A) is 
commutative, then efficient determinant computation is feasible and otherwise determinant is 
as hard as permanent. 
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1 Introduction 



Given a matrix M = {mjj}, the determinant of M, denoted by det(M) is given by the polynomial 
det(M) = X]o-es„ sgn(cr) fl"^^ m-jo-j, while the permanent of M, denoted by per(M) is defined by 
the polynomial per(M) = X]o-e5,i HiLi '^^o"*- Though deceivingly similar in their definitions, the 
determinant and permanent behave very differently with respect to how efficiently one can compute 
these quantities. The determinant of a matrix over any field can be efficiently computed using 
Gaussian elimination. In fact, determinant continues to be easy even when the entries come from 
some commutative algebra, not necessarily a field [Sam42, Ber84, Chi85, MV97]. Computing the 
permanent of a matrix over the rationals, on the other hand, as famously shown by Valiant [Val79], 
is just as hard as counting the number of satisfying assignments to a Boolean formula or equivalently 
T^P-complete even when the entries are just and 1. Given this state of affairs, it is natural to ask, 
"what is it that makes the permanent hard while the determinant is easy?" Understanding this 
distinction in complexity of computing the determinant and permanent of a matrix is a fundamental 
problem in theoretical computer science. 

Nisan first pioneered the study of noncommutative lower bounds in his 1991 groundbreaking 
paper [Nis91]. In one of that paper's more important results, Nisan proves that any algebraic 
branching program (ABP) that computes the determinant of a matrix M = {rriij} over the non- 
commutative free algebra F(a;ii,. . . ,Xnn) niust have exponential size; this then implies a similar 
lower bound for arithmetic formulas. This contrasts markedly with the many known efficient algo- 
rithms for determinant in commutative settings, which include polynomial-sized ABPs [MV97]. 

This problem takes on added significance in light of a connection discovered by Godsil and 
Gutman [GG81] and developed by Karmarkar et al. [KKL"*"93] between computing determinants 
and exponential time algorithms for approximating the permanent. The promise of this approach 
was cemented when Chien et al. [CRS03], expanding on work by Barvinok [Bar99], showed that if 
one can efficiently compute determinant of an n x n matrix M whose entries rriij are themselves 
matrices of O(n^) dimension, then there is a fully polynomial randomized approximation scheme 
for the permanent of a 0-1 matrix; similar results were later proven by Moore and Russell [MR09]. 
Thus understanding the complexity of noncommutative determinant is of both algorithmic and 
complexity-theoretic importance. 

Nisan's results are somewhat limited in that they apply only to the free algebra F(xj) and not 
to specific finite dimensional algebras (such as those used to approximate the permanent), and 
because they do not apply outside of ABPs and arithmetic formulas. Addressing the first concern, 
Chien and Sinclair [CS07] significantly strengthened Nisan's original lower bounds to apply to a 
wide range of other algebras by analyzing those algebras' polynomial identities. In particular, they 
show that Nisan's lower bound extends to d x d upper-triangular matrix algebra over a field of 
characteristic for any d > 1 (and hence over Mrf(F), the full d x d matrix algebra as well), the 
quaternion algebra, and several others, albeit only for ABPs. 

In a significant advance, Arvind and Srinivasan [ASIO] recently broke the ABP barrier and 
showed noncommutative determinant lower bounds for much stronger models of computation. 
They show that unless there exist small circuits to compute the permanent, there cannot exist 
small noncommutative circuits for the noncommutative determinant. More devastatingly from the 
algorithmic point of view, they show that computing det(M) where the rriij are linear-sized matrix 
algebras is at least as hard as (exactly) computing the permanent. Arvind and Srinivasan thus bring 
into serious doubt whether the determinant-based approaches to approximating the permanent are 
computationally feasible. 
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While these collections of results make substantial progress in our understanding of when de- 
terminant can be computed over a noncommutative algebra, they are still incomplete in significant 
ways. First, we do not know whether Arvind and Srinivasan's results rule out algorithms for de- 
terminants over constant-dimensional matrix algebras, which are still of use in approximating the 
permanent. More expansively, we still do not know the answer to what is perhaps the fundamental 
philosophical question underlying this: 

Whether there is any noncommutative algebra over which we can compute determinants 
efficiently, or whether, as may seem attractive, commutativity is a necessary condition 
to having such algorithms? 

1.1 Our results 

In this paper, we fill in most of these remaining gaps. Our first main result extends Arvind and 
Srinivasan's results all the way down to 2 x 2 matrix algebras. 

Theorem 1.1. (stated informally^ ) Let M2(F) be the algebra of 2x2 matrices over a field F. Then 
computing the determinant over M2{¥) is as hard as computing the permanent over F. 

The proof of this theorem works by retooling Valiant's original reduction from #3SAT to per- 
manent. One would not expect to be able to modify Valiant's reduction to go from ^^SSAT to 
determinant over a field F, as there are known polynomial-time algorithms in that setting. How- 
ever, when working with M2 (F) , what we show is that there is just enough noncommutative behavior 
in M2(F) to make Valiant's reduction (or a slight modification of it) go through. 

Given the central role of matrix algebras in ring theory, this allows us to prove similar results for 
other large classes of algebras. In particular, consider a finite-dimensional algebra A over a finite 
field F. This algebra has a radical R{A), which happens to be a nilpotent ideal of A. Combined 
with classical results from algebra (in particular the simple structure of the 2x2 matrix algebras) 
the above theorem can be extended as follows to yield our second main result. 

Theorem 1.2. (stated informally^) If A is a fixeS finite dimensional algebra over a finite field 
such that the quotient A/R[A) is noncommutative, then computing determinant over A is as hard 
as computing the permanent. 

In particular, if the algebra is semisimple (i.e, R{A) = 0), then the commutativity of A itself 
is determinative: if A is commutative, there is an efficient algorithm for computing det over A] 
otherwise, it is at least as hard as computing the permanent. The class of semisimple algebras 
includes several well-known examples, such as group algebras. 

It may be tempting at this point to see the sequence of lower bounds starting from Nisan's 
original work and conjecture that computing det over A for some algebra A is feasible if and only 
if A is commutative. Perhaps surprisingly, we show that this is not the case — in, fact there do 
exist noncommutative algebras A for which there are polynomial-time algorithms for computing 
det over A. For instance, in our third main result, we show that computing the determinant where 

^See Theorem 3.5 for a formal statement. 
"'See Theorem 5.1 for a formal statement. 

''By fixed, we mean that the algebra is not part of the input; we fix an algebra A and consider the problem of 
computing the determinant over A. 
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the matrix entries are d x d upper triangular matrices for constant d is easy. For reasons that will 
soon be clear, we will state this result, more generally, in the language of radicals. 

Theorem 1.3. Given a finite dimensional algebra A and its radical R{A), let d be the smallest 
value for which R{A)'^ = (i.e. any product of d elements of R{A) isO). If A/R[A) is commutative, 
there is an algorithm for computing det over A in time poly(n'^). 

While this description of the class of algebras that allow efficient determinant computation is 
somewhat abstruse, it does include several familiar algebras. Perhaps most familiar is the algebra 
C/d(F) of d X d upper-triangular matrices, for which R{Ud{¥))'^ = 0. What the result states is 
that the key to whether determinant is computationally feasible is not commutativity alone. For 
noncommutative algebras, it is still possible that determinant can be efficiently computed, so long 
as all of the noncommutative elements belong to a nilpotent ideal and have a limited "lifespan" of 
sorts. 

The above theorems together yield a nice dichotomy for constant dimensional algebras over a 
finite field. Given any such algebra A of constant dimension D over a finite field, either A/R{A) is 
commutative or not. Furthermore, if A/R{A) is commutative, we have that R{A) is nilpotent with 
nilpotency index at most D which is a constant. We thus, have the following dichotomy: if A/R{A) 
is commutative, then efficient determinant is feasible else determinant is as hard as permanent. 

Does this yield a complete characterization of algebras over which efficient determination com- 
putation is feasible? Unfortunately not. In particular, what if the dimension D is non-constant, i.e., 
the algebra is not fixed but given as part of the input or if the algebra is over a field of characteristic 
0? In these cases, the lower bound of Theorem 1.2 and upper bound of Theorem 1.3 are arguably 
close, but do not match. A complete characterization remains an intriguing open problem. 

Organization of the paper: After some preliminaries in Section 2, we prove lower and upper 
bounds in two concrete settings: we prove a lower bound for 2x2 matrix algebras in Section 3 and 
an upper bound for small-dimensional upper triangular matrix algebras in Section 4. The results 
on general algebras are in Section 5, followed by some discussion in Section 6. 

2 Preliminaries 

In this section we define terms and notation that will be useful later. 

An (associative) algebra A over a field F is a vector space over F with a bilinear, associative 
multiplication operator that distributes over addition. That is, we have a map ■ : A x A ^ A that 
satisfies: (a) x ■ {y ■ z) = (x • y) • z for any x,y,z £ A, (b) A(x • y) = (Ax) - y = x ■ (Ay), for any A € F 
and x, ?/ G j4, and {c) x ■ {y + z) = x ■ y + x ■ z and {y + z) ■ x = y ■ x + z ■ x ior any x,y, z £ A. We 
will assume that all our algebras are unital, i.e., they contain an identity element. We will denote 
this element as 1. For more about algebras, see Curtis and Reiner's book [CR62]. A tremendous 
range of familiar objects are algebras; we will be concerned with the algebra of dx d matrices over 
F, which we will denote M(i(F), as well as the algebra of d x d upper-triangular matrices over F, 
or C/rf(F). Other prominent examples are the free algebra F(xi), the algebra of polynomials F[xj], 
group algebras over a field, or a field considered as an algebra over itself. 

Given an n x n matrix M = {rriij) whose elements belong to an algebra A, the determinant of 
M, or det(M), is defined as the polynomial det(M) = X^o-es^ sgn((T) mjo-j. Note that when 
A is noncommutative, the order of the multiplication becomes important. When the order is by 
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row, as above, we are working with the Cayley determinant. The permanent of the same matrix 
is per(M) = J2aeSn TYi=i ^i^i- We will denote by det^i (and per^^) the problem of computing the 
determinant (and permanent) over an algebra A. 

We recall also the familiar recasting of the determinant and permanent in terms of cycle covers 
on a graph. Suppose M = (niij) is an n x n matrix over an algebra A. Let G{M) denote the 
weighted directed graph on vertices 1, . . . ,n that has M as its adjacency matrix. A permutation 
vr : [n] — 7- \n\ from the rows to the columns of M can be identified with the set of edges (z,7r(i)) 
in the graph G{M); it is easily observed that these edges form a (directed) cycle cover of G{M). 
Letting C{G) denote the collection of all cycle covers of G{M), we can write 

det(M) = J2 sgn(C)mi^c'(i)"i2,c(2) • • -"in.cH (^-l) 

CeCiG{M)) 

and 

per(M) = ^ mi^c{i)m2,c{2)---mn,c{n), (2-2) 

CeC(G{M)) 

where for a given cycle cover C, C{i) represents the successor of vertex i in C, and sgn(C) is the 
sign of G. It is known that sgn(C) = (—1)""'^, with c being the number of cycles in C, and that this 
is also the sign of the corresponding permutation. We will denote the weight of an edge e = {x, y) 
as w{e) or w{x, y). Further, for a subset of edges B = {{xi,yi), . . . , (x|5|, y\B\)} of a cycle cover C 

\B\ 

with Xi < Xj+i, we can define the weight of B as w{B) = ni=i w{xi,yi). (Note that the product is 
in order by source vertex.) Thus w{C) = J|- rrii^cii) by w{G) is the weight of the cycle cover, and 
the product sgn(C) fnifi(i) is the signed weight of G. 

3 The lower bound for 2x2 matrix algebras 

In this section, we show our key lower bound for 2x2 matrix algebras. Our proof is based on 
Valiant's seminal reduction from #3SAT to permanent, as modified by Papadimitriou [Pap94] and 
also described in the complexity textbook by Arora and Barak [AB09]. We first give a self-contained 
description of that, before detailing our modifications of it. 

3.1 Valiant's lower bound for the permanent 

Valiant's reduction is from #3SAT to permanent; given a #3SAT formula (p onn variables and m 
clauses, he constructs a weighted directed graph G^p on poly(n, ra) vertices such that the number 
of satisfying assignments of is equal to constant x per(M(G<^)), where M{G^p) is the adjacency 
matrix of Gip. The key components of are the variable, clause, and XOR gadgets shown in 
Figure 1.^ The idea is that there will be a relation between satisfying assignments of Lp and cycle 
covers of G,p\ moreover, for each satisfying assignment, the total weight of its corresponding cycle 
covers will be the same. 

Before defining G^p itself, we first work with a preliminary graph that contains n variable 
gadgets and m clause gadgets, but no XOR gadgets; all of the gadgets are disjoint from each other. 

^We follow a convention from [AB09] in allowing gadgets to sometimes have multiple edges between the same two 
vertices. While technically prohibited in a graph defined by a matrix, this can be fixed by adding an extra node in 
these edges. 



4 



-1 




variable clause XOR 



Figure 1: Gadgets used in proving Valiant's lower bound. All edges have weight 1 unless noted otherwise. 
For the variable and clause gadgets, the solid (not dotted) edges are called (vertex) external edges or (clause) 
external edges. Note that the number of external edges in a variable gadget is not fixed, and need not be 
the same for the True and False halves of the gadget. 



For the moment, the number of external edges in each of the variable gadgets is unimportant. In 
analyzing G^, we will use the following: 

Lemma 3.1. The following hold for the gadgets in Figure 1: (a) A variable gadget has exactly two 
cycle covers. Each cycle cover contains one long cycle using all of the external edges on one side 
of the gadget and the long middle edge, as well as all the self-loops on the other side of the gadget, 
(h) In a clause gadget, there is no cycle cover that uses all three external edges. For every proper 
subset S of the external edges in a clause gadget, there is exactly one cycle cover that contains 
exactly the edges in S; this cycle cover has weight 1. 

As all n + m gadgets in are disjoint, any cycle cover of will be a union of n + m smaller 
cycle covers-namely, one for each gadget. The choice of cycle cover for each gadget defines the 
value of each variable and which literals are satisfied in each clause. 

More precisely, for a variable gadget, let the term True cycle cover denote the cycle cover 
containing the external edges on the True side of the gadget. Analogously, the False cycle cover 
refers to the cycle cover containing the external edges on the False side of the gadget. The idea 
is that a cycle cover of G^ sets a variable to T or F by choosing either the True or False cycle 
cover. Meanwhile, for clause gadgets, the intention is that each external edge will correspond to 
one of the three literals in the clause, and an external edge is used in a cycle cover if and only if 
the corresponding literal is set to F (i.e. the corresponding literal is not satisfied). Since no cycle 
cover can contain all three external edges of a clause gadget, in this interpretation at least one of 
the literals in the clause must be satisfied. 

We say a cycle cover C of is consistent if (1) whenever C contains the True cycle cover of 
the gadget for a variable x^, it contains all clause external edges for instances of the negative literal 

and no clause external edges for instances of the positive literal Xk, and (2) conversely, whenever 
C contains the False cycle cover for Xk, it contains all clause external edges for instances of x^ but 
no clause external edges for instances of xj^. A consistent cycle cover therefore does not "cheat" by 
claiming to set Xk to T (for example) in a variable gadget but to F in a clause gadget. This is close 
to what we want: 
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Lemma 3.2. The number of satisfying assignments of <p is equal to the total weight of consistent 
cycle covers of . 

Proof. This follows from combining the natural bijection between satisfying assignments and con- 
sistent cycle covers and the fact from Lemma 3.1 that every cycle cover of a clause gadget has 
weight 1. □ 

Of course, nothing about guarantees that a cycle cover must be consistent, and in fact many 
inconsistent covers exist. To fix this, we need to use the critical XOR gadgets to obtain the final 
graph G^. 

The graph G^ is constructed as shown in Figure 2 (left). It has the same n variable gadgets and 
m clause gadgets as G^, with the gadget for each variable having as many True external edges 
as there are instances of Xk in (/?, and as many False external edges as there are instances of lEk. 
Now, however, for each appearance of a literal Xk or x^" in a given clause, an XOR gadget is used 
to replace the corresponding external edge in that clause gadget and a distinct external edge on 
the appropriate side of the variable gadget for Xk- The role of the XOR gadgets is to neutralize the 
inconsistent cycle covers of while still maintaining the property that each satisfying assignment 
of if contributes the same to the total weight of cycle covers. This leads to the description of the 
final graph G^ itself. 

We now state the important properties of the XOR gadget, the key component of Valiant's 
proof. 

Lemma 3.3. Suppose a graph G contains edges {u,u') and (v,v'), with all four vertices distinct. 
Suppose now that the edges {u,u') and {v,v') are replaced by an XOR gadget as shown in Figure 1, 
resulting in a new graph G' (with four new vertices a, b, c and d). Let C„y^ be the set of cycle covers 
containing {u,u') but not {v,v'), and w^^^ = X^ceC, \ ^(C*) be their total weight. Let C„\^ and 
Wy = X^c'gc ^ ^(C) ^6 defined analogously. Then there exist two disjoint sets of cycle covers of G' 
with total weight 4:Wu\v o-nd 4tf„\„, while all cycle covers of G' not in these sets have total weight 
0. 

The proof is omitted, as we will state and prove our own modified version of this in Section 3.2. 
This leads to the following: 

Theorem 3.4. [Valiant] Given a 3-SAT formula (p and the graph G^ as described, per(G,^) = 
4'^™ 5", where S is the number of satisfying assignments of ip. 

We omit the formal proof, but give some of the intuition. Beginning with G^, we begin adding 
XOR gadgets one at a time. When a pair of edges is replaced by an XOR gadget, any cycle covers 
that are consistent with respect to that pair of edges are turned into a set of cycle covers whose 
total weight is a factor of 4 more than the original weight. All other cycle covers in the new graph 
have total weight 0. This continues until each of the 3m XOR gadgets are added, at which point 
the original consistent cycle covers have become a set of cycle covers with total weight 4^™ while 
all other cycle covers in the final graph have weight 0. The total weight of the cycle covers in the 
final graph is therefore 4^™'5, as required. 

3.2 Our construction 

We now prove the following: 
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Figure 2: Left: Subgraph of corresponding to clause (a;i V X2 V xl), with the clause gadget in center. 
Three variable gadgets are connected to the clause gadget via XOR gadgets. Right: Examples of how 
gadgets may have cycle covers of different sign. 



Theorem 3.5. Let ¥ he a field of characteristic p > 0. If p = 0, computing det^,/2(F) i^P-hard. 
On the other hand, if p > and odd, then computing detjv/2(F) 'is ModpP-hard. 

Our proof is also a reduction from #3SAT (or Modp-SAT in the case of positive odd character- 
istic) and is based on Valiant's framework as described in the previous subsection. Given a 3SAT 
formula ip, we wish to construct a directed graph H^p with weights belonging to M2(F) such that 
the number of satisfying assignments of (p can be computed from det{M{H^)), as expressed in 
equation (2.1) above. We will first describe the graph and then prove its correctness. 

A very naive but instructive first try would be to simply use the graph from Valiant's 
construction, replacing each edge weight w £ ¥ with wl, where I2 is the 2x2 identity matrix. 
This fails, of course, because of the sign factor sgn((7) inside the summation, which is based on 
the parity of the number of cycles in C. The immediate problem is that each of the three types 
of gadgets could conceivably use an odd or even number of cycles. As shown in Figure 2 (right), 
variable gadgets may have a different number of self-loops on different sides; clause gadgets may 
use one or two cycles depending on which external edges are chosen; and XOR gadgets show similar 
behavior. 

Fortunately, these problems can be overcome if we also allow ourselves to modify the edge 
weights, and crucially, use the noncommutative structure available in M2(F). This results in the 
gadgets shown in Figure 3. We now define two graphs, a preliminary graph and final graph 
H^, in analogy with and from Section 3.1. The new graphs and H^^ will be constructed 
in the same manner as G^p, only using the modified gadgets from Figure 3 instead of the original 
gadgets in Figure 1. 

The rough idea behind these gadgets is that with the new weights, each resulting cycle cover of 
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a gadget of the "wrong" sign will have an extra —1 sign from its edge weights. The determinant is 
then essentially the same as the permanent. We now explain the changes in more detail. 

For variable gadgets, the fix is easy - all we have to do is make sure that both sides of the 
gadget have an even (for example) number of vertices, and hence an even number of self loops. 
This can be accomplished by adding, if necessary, a new vertex and appropriate new edges on one 
or both sides. The new external edges, if any, will not be connected to any of the clause gadgets. 

For clause gadgets, we need to address the problem that some cycle covers have only one cycle, 
while others have two. Here we benefit from the observation that one of the edges, (x, y) in Figure 3, 
is used only in cycle covers with two cycles. Thus we can correct for parity by changing the sign 
of this edge from I2 to —I2; as a result, every cycle cover of a clause gadget has the same signed 
weight. 

For XOR gadgets, simply changing the edge weights to scalar multiples of I2 is insufficient. 
(Indeed, Valiant presciently noticed this in 1979!) However, we can save the construction by using 
more sophisticated matrix-valued edge weights instead. In particular, we define the following three 
2x2 matrices: 

--{I -"0^^-K-> 'oy^-n -:)■ 

We then modify the weights of the edges between vertices a, b, c and d. Specifically, each edge 
entering vertex b has its weight multiplied by X; each edge entering c has its weight multiplied by 
Y, and each edge entering d has its weight multiplied by Z. 

For now, with i/^ defined, we prove that computing det{H^) is equivalent to computing the 
number of satisfying assignments of (p. We first observe the following analogue of Lemma 3.2. 

Lemma 3.6. Let Ccon be the set of all consistent cycle covers of H^. Then there exists z G {1,-1} 
such that for all C G Ccon; we. have sg'n.{C)w{C) = zl2- 

Proof. As in the proof of Lemma 3.2, there is a bijection between satisfying assignments of if and 
consistent cycle covers of H^. We need to show that each of these cycle covers has the same signed 
weight. For such a cycle cover C S Ccon we have sgn(C) = (— 1)"'h~'^('^), where is the number of 
vertices in and c(C) is the number of cycles in C. We further know that (— 1)^^'-') = (— 1)?+'^+'?^ 
where p is the number of cycles used to cover the n variable gadgets, m is the number of clauses. 
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and q is the number of times C uses two cycles to cover a clause gadget. Since we assumed p to be 
even, we have sgn(C) = (— 1)"-//+™-+'?. 

On the other hand, w{C) is the product of the edge weights of C. All of these weights are 
I2 except for the w{x,y) in the clause gadget, which has weight —I2 and shows up when C uses 
two edges for a clause gadget. Thus w{C) = {—1)^I2, and sgn{C)'w{C) = (—!)"« +™/2) which is 
independent of the cycle cover C. (Hence, J^ceCcon^^^C^^^i^^ = (-l)"^+™'5/2, where S is the 
number of satisfying assignments of 93). □ 

Without loss of generality, we can assume from here on that the sign z is positive, as we can 
insert a new vertex within an edge so that + m is even. 

We now prove the following useful identities of XOR gadgets, which can be verified by hand: 

Lemma 3.7. Let MxoR be the adjacency matrix for the XOR gadget, or 



Letting Mij indicate the minor of M with row i and column j removed, we have (1) det(M3^i) = 
-4/2, (2) 'det(Mi,3) = -4J2, (3) det(M) = det(Mi,i) = det(M3,3) = det(Mi3,i3) = 0, where 



Now consider a graph G with vertices labeled l,...,nG and weights in M2(F). Suppose G 
contains vertex-disjoint edges {u,u') and {v,v'), each with weight 12- Suppose now that the edges 
{u,u') and {v,v') are replaced by an XOR gadget as shown in Figure 3. This results in a new graph 
G' , with four new vertices a, b, c and d, which we number hq + 1, . . . , no + 4. We now define a 
mapping tp from C{G) to subsets of C'(G) as follows: Given cycle covers C G C(G) and C' £ C{G'), 
then C' £ ii){C) if and only if (1) for all edges e G C\{(n, u'), (u, v')}, we have e G C", (2) (-u, u') G G 
if and only if (n, a), (c, u') G C", and (3) {v, v') G C if and only if (f , c), (a, v') G C". 

This leads to the following analogue of Lemma 3.3: 

Lemma 3.8. Let Cu\v = {C G C(G) : {u,u') G G,{v,v') C} be the set of cycle covers of 
G containing {u,u') but not {v,v'), and Cy\u = {C G C{G) : {v,v') G C,{u,u') C}. Then there 
exists a mapping Tp from C{G) to subsets ofC'{G) such that ip{Ci)riip{G2) = for all Ci, C2 G C{G) 
and (1) for any C G Cu\v, the total weight ofip{C) is ^^(^/^^(c) sgn(C")i(7(C) = 4sgn{G)w{C), (2) 
for any G G C^\u, Ylc'&ip{C) sgn(C)w(C") = 4sgn(C)w(C) J2, and (3) the remaining cycle covers in 
G' have total weight Ec'^^(C)VC6C„\„uc„\„ sgn{G')w{G') = 0. 

Proof. We start with proving (1). Fix any G G Cu\^. Notice that 4'{G) consists of all G' G C{G') 
that contain {u,a), {c,u') and all of C's edges except (u,u'). Call this set of common edges Ec; 
by the assumption that w{u,u') = L2, w{Ec) = w{G) The set ipiG) consists of all possible ways of 
completing Eq to a cycle cover C" of G' by adding edges to G' so that every vertex has indegree 
and outdegree 1. Within Ec, the only vertices with deficient degree are a,b,c and d. Vertices b 
and d have indegree and outdegree 0, while a has indegree 1 and outdegree 0, and c has indegree 
and outdegree 1. Note that the edges (n, a) and {c,u') must belong to the same cycle in G' , and 
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so the edges in Ec form zero or more completed cycles and an incomplete cycle from c to a. The 
nmnber of completed cycles is c(C) — 1, where c(C) is the number of cycles in C. 

We thus need to add three edges matching the vertices {a, b, d} to the vertices {b, c, d}; call these 
three edges -BxoR) so that E'cUE'xoR forms a cycle cover C . The weight of C is therefore w{C') = 
w{Ec)w{Exor)- The sign of C is (— where c(C") is the number of cycles in C . We can 
see that c(C") is the sum of the number of completed cycles in Eq and the number of cycles among 
{a, b, c, d} assuming the existence of an edge from c to a. Hence c(C") = c(C) — 1+c(SxorU{(c, a)}, 
and so sgn(C7') = -sgn(C) (-1)^-'=(-^xorU{ (c,a)}) ^ _sgn(C)sgn(ExoR U {(c, a)}). 

Thus, Ec'e^{C) sgn{C')w{C') = -sgn{C)w{Ec) J2c'&ip(C) sgn(£;xoR U {(c, a)})u;(^xoR) = 
— sgn{C)'w{Ec) det(M3^i). From Lemma 3.7, this is 4sgn{C)'w{Ec) = ^sgn{C)w{C), as required. 

The proof of (2) proceeds similarly, except that ipiC) contains {v, c) and (a, v') instead of {u, a) 
and (c, v'). The set of common edges then has an incomplete path from a to c. As a result, we end 
up with Ecey^iO ^MC')w{C') = -sgn{C)w{Ec)det{Mi,3) = Asgn{C)w{C)J2. 

To prove (3), we observe that a cycle cover in C{G') that contains (n, a) and (c, u') but not {v, c) 
or (a, v') must fall into Tp{C) for some C € C„\„; similarly, any cycle cover containing {v, c) and (a, u') 
but not {u, a) or c, u') must fall into ip{C) for some C G C„\„. These were already accounted for in 
the proofs of (1) and (2), so we can concentrate only on the leftover cycle covers. Partition these 
leftover cycle covers into equivalence classes based on their edge sets excluding edges wholly within 
{a, 6, c, d}; namely C[ ~ C2 if and only if C( \ {a, b, c, d} x {a, b, c, d} = C2 \ {a, b, c, d} x {a, 5, c, d}. 
For any equivalence class, its cycle covers must either all (a) contain none of these four edges, (b) 
contain (u, a) and (a, v') only, (c) contain {v, c) and (a, v') only, or (d) contain all four edges. 

Up to sign, the total weights of those equivalence classes in (a) contain a factor of det(M), those 
in (b) contain a factor det(Mi^i), those in (c) contain a factor det(M3^3), and those in (d) contain a 
factor det(Mi3_i3). From Lemma 3.7, all four of these determinants are 0, and so the total weights 
of the cycle covers in any equivalence class is 0, as is therefore the total weight of all the leftover 
cycle covers. □ 

With this in hand, we can prove the key result: 

Theorem 3.9. Given a 3SAT formula ip with S satisfying assignments, let the graph with 
weights in M2(F) be as defined above. Then det{Hi^) = a/2 + 6J2, where a + b = 4^™5. 

Proof. The structure of the proof is similar to the sketch given after Theorem 3.4, though with 
extra care needed for the complications of working with matrices. In the end, each cycle cover of 
G ends up with weight 4'^™'/2 or 4'^'" J2, giving the result. 

Let us start with H^, which we know from Lemma 3.6 has det{H^) = Sl2. In particular, for 
each satisfying assignment of 93, there is a consistent cycle cover of i/^ of weight l2- There exist 3m 
pairs of edges in that when replaced by XOR gadgets will convert to H^; each of these pairs 
contains an external edge in a clause gadget and an external edge in a variable gadget referring to 
the same literal. 

Consider what happens when we replace one of the above-mentioned edge pairs with an XOR 
gadget, forming a new graph H^. From Lemma 3.8, each cycle cover C that is consistent on 
this edge pair in will be mapped to tp{G), a set of cycle covers in the new graph whose total 
signed weight will either be 4/2 or 4J2. Further, since all of these sets ip{G) are disjoint and all 
other cycle covers have total signed weight 0, the total signed weight of all cycle covers in is 
J2ceC^ (G) 4/^2 (C*), where C^^^{G) are those cycle covers of G that are consistent on this edge pair, 
and K2{G) is either I2 or J2. 
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Now suppose a second edge pair is replaced with an XOR gadget, resulting in the graph H^. 
Consider a cycle cover C of -ff^ in ^(C) that is consistent on both the first and second edge pairs. 
Then each cycle cover of in ipiC) will be mapped to a set of cycle covers il){il){C)) of i?^, with 
signed weight that is 4/2 or 4J2 multiple of its signed weight in H^. The set ip{iJj{C)) therefore 
has total signed weight of either I6/2 or I6J2, since all of the images of ip are disjoint. Once again, 
the total signed weight of all cycle covers in is 'YliceC^''^ {G) 16-^2(C)i where Ccon{G) is the set of 
cycle covers of G consistent on both edge pairs. 

Carrying this out over all 3m edge pairs to reach i?^, we see that every consistent cycle cover 
of becomes a disjoint set of cycle covers in H^p of total signed weight 4^'^/2 or A^"^J2, while all 
other cycle covers in have total weight 0. The total weight over all original consistent cycle 
covers is X^ceCconCC) 4:'^"^K2{C). This therefore takes the form given in the theorem. □ 

This completes the proof of Theorem 3.5. 

4 Computing the determinant over upper triangular matrix algebras 

In this section, we consider the problem of computing the determinant over the algebra of upper 
triangular matrices of dimension d. We show that the determinant over these algebras can be 
computed in time N'^^'^\ where N denotes the size of the input. We will then generalize this 
theorem to arbitrary algebras to yield Theorem 5.5. 

Given a field F, we denote by f/d(F) the algebra of upper triangular matrices of dimension d 
with entries from the field F. 

Theorem 4.1. Let¥ be afield. There exists a deterministic algorithm, which when given as input 
an n X n matrix M with entries from [/^(F), computes the determinant of M in time poly(A^'*), 
where N is the size of the input. 

Proof. The algorithm is simple. We write out the expression for the determinant of M and note 
that each entry of det(M) may be written as the sum of n'^^'^^ many determinants of matrices with 
entries from the underlying field. Since each of these can be computed in time N'-'^^\ we obtain 
an A^'^('^)-time algorithm for our problem. 

Let M = {rnij)ij, where ruij £ Ud{¥) for each i,j £ [n]. Given m £ Ud{¥), we use m(p,q) to 
denote the {p, q)th. entry of m. We have 

det(M) = ^ sgn(a)mi_^(i)m2,^(2) • • • m„_^(„) 

Consider a product of matrices m = mi ■ ■ ■ rUn where each mi G [/^(F). For p,q £ [d] such that 
p < q, we may write the {p, q)th entry of m as 

m{p,q)= ^ mi{p,ki)m,2{ki,k2) ■ ■ ■mn{kn-i,q) 

ki,k2,...,kn-i£[d] 

= ^ mi{p,ki)m2{ki,k2) ■ ■ ■mn{kn-i,q) (4.1) 

p<ki<---<kn-i<q 

where the last equality follows since mi{k,l) = unless k < I. Note that the number of terms in 
the summation in (4.1) is equal to the number of increasing sequences of length n consisting of 
elements from [d] and is bounded by n^^'^\ 
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Fix any p,q £ [d] such that p < q. By (4.1), we may write det(M)(p, q) as 
det{M){p,q) = ^ ^ sgn(cr) • /ci) • m2,^(2)(A;i, ^2) • • • m„^^(„)(/c„_i, g) (4.2) 

p<ki<---<k„-i<q cr&Sn 

We now note that each of the inner summations may be written as the determinant of an 
appropriate matrix over the underlying field. Fix any k = (ki, . . . , /c„_i) satisfying p < ki < k2 < 
■ ■ ■ < kn^i < q. Denote by the matrix {mij{ki-.i, ki))ij, where ko denotes p and ki denotes q. 
It follows from (4.2) that det{M){p,q) = Efcdet(A%). 

Note that the matrices are nx n matrices with entries from the underlying field and hence, 
their determinants can be computed in time N^^^\ Hence, we can compute det{M)(p,q) — for 
each p,q — in time n'^('^) • iVO(i) = NOid) _ -pj^g result follows. □ 

5 Determinant computation over general algebras 

We now consider the problem of computing the determinant of an n x n matrix with entries 
from a general finite-dimensional algebra A of dimension D over a field F that is either finite, or 
the rationals. We consider two algorithmic questions: the first is the problem of computing the 
determinant over A, where A is a fixed algebra (and hence of constant dimension) such as M2(F); 
the second is the case when A is presented to the algorithm along with the input (in this case, A 
could have large dimension). We present our results for the latter case in the appendix. 

In the first case, we prove a strong dichotomy for finite fields of characteristic p > 2. For any 
fixed algebra A, we show, based on the structure of the algebra, that either the determinant over 
A is polynomial-time computable, or computing the determinant over A is ModpP-hard. 

We first recall a few basic facts about the structure of finite dimensional algebras. An algebra 
is simple if it is isomorphic to a matrix algebra (possibly of dimension 1) over a field extension of 
F. An algebra is said to be semisimple if it can be written as the direct sum of simple algebras. " 

Recall that a left ideal in an algebra yl is a subalgebra / of ^ such that for any x G I and a G A, 
we have ax € I; a right ideal is defined similarly. An ideal / is said to be nilpotent if there exists 
an m > 1 such that the product of any m elements from / is 0. The radical of A denoted R{A) 
is defined to be the ideal generated by all the nilpotent left ideals of A. We list some well-known 
properties of the radical (see [CR62, Chapter IV]): (a) The radical is a left and right ideal in A, 
(b) The radical is nilpotent: that is, there exists a d G N such that the product of any d elements 
of R{A) is 0. The least such d is called the nilpotency index of the radical R{A), and (c) A/R{A) 
is semisimple. 

An algebra A is a semidirect sum of subalgebras Bi and i?2 if ^ = -Bi © B2 as a vector space; 
we denote this as A = Bi ®'i?2- The Wedderburn-Malcev theorem (Theorem A. 3) tells us that any 
algebra is a semidirect sum of its radical with a subalgebra. We refer to such a decomposition as a 
Wedderburn-Malcev decomposition. 

We start with the hardness result. 

Theorem 5.1. Let A denote any fixed algebra over a finite field F of characteristic p > 2. If 
A/R{A) is non-commutative, computing the determinant over A is ModpP-hard. 

"This is not the standard definition of semisimphcity in the case of infinite fields. However, we will only use it in 
the case that F is finite. See Appendix A. 
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Proof. Consider the problem of computing the determinant over an algebra A such that A/R{A), 
the "semisimple part" of A, is non-commutative. Since A/R{A) is semisimple, we know that 
A/R{A) = where each Ai is a simple algebra, and hence isomorphic to a matrix algebra 

over a field extension of ¥. If each of the AiS is a matrix algebra of dimension 1 (that is, each Ai 
is simply a field extension of F), then A/R{A) is commutative. Hence, w.l.o.g., we assume that Ai 
has dimension greater than 1. Moreover, by the Wedderburn-Malcev theorem (see Theorem A. 3 
in the appendix), we know that A contains a subalgebra B = A/R{A). Thus, the algebra Ai 
is isomorphic to a subalgebra of A. Thus, Theorem 3.5 immediately implies that computing the 
determinant over A is ModpP-hard. □ 

5.1 The upper bound 

In this section, we show that if A/R{A) is commutative, then the determinant over A is efficiently 
computable. However, we present our result in some generality, which will be useful later. We 
assume that the algebra A is presented to the algorithm along with the input as follows: we are 
given a (vector space) basis {ai, . . . , ud} for A along with the pairwise products aiUj for every 
[D]. Let d denote the nilpotency index of R{A). 

The Wedderburn-Malcev theorem (see Theorem A. 3) tells us that the algebra A = B (B' R{A), 
where S is a semisimple subalgebra of A isomorphic to A/R{A), and hence commutative. 

We will use without explicit mention the following result, which was explicit in the work of 
Chien and Sinclair [CS07], and implicit in that of Mahajan and Vinay [MV97] (and also many 
other works): 

Theorem 5.2. There is a deterministic algorithm which, when given any commutative algebra A 
of dimension D and an n x n matrix over A as input, computes the determinant of A in time 
poly{n,D). 

We start with two simple lemmas. 

Lemma 5.3. There is a deterministic polynomial-time algorithm which, when given an algebra A 
as input, computes the nilpotency index of A. 

Proof. Let d denote the nilpotency index of A. It is easy to see that d < D, the dimension of 
the algebra as a vector space over F. The algorithm computes a basis for R{A) (this can be done 
in deterministic polynomial time by Theorem A.l), and then successively computes a basis for 
R{A)'^, R{Af, R{A)^ and outputs the least d such that R{A)'^ = {0}. □ 

Lemma 5.4. Let A be a finite- dimensional algebra with Wedderburn-Malcev decomposition A = 
B®' R{A). Then, le B. 

Proof. We can write the identity 1 of j4 as 1 = 6 + r, where b £ B and r E R(A). We would like to 
show that r = 0. Note that b = b ■ 1 = b^ + br. Since £ B and br £ R, we must have br = 0. 
Similarly, r = Ir = 6r + r^. But br = implies that r = r^. This implies that r = r^ for any k > 1. 
But we know that r is nilpotent. Hence, r = 0. □ 

These lemmata and a generalization of Theorem 4.1 yield the following: 
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Theorem 5.5. There exists a deterministic algorithm, which when given as input an algebra A of 
dimension D s.t. A/R{A) is commutative and an n x n matrix M with entries from A, computes 
the determinant of M in time N'^'^'^\ where d is the nilpotency index of R{A) and N is the size of 
the input. 

In particular, when A is a fixed algebra, then d < D = 0(1), and hence, Theorem 5.5 gives us 
a polynomial-time algorithm. This yields straightaway the sharp dichotomy theorem in the case of 
a fixed algebra over finite fields of odd characteristic. 

Corollary 5.6. Let F he any finite field of odd characteristic and A he any fixed algebra over 
F. Then, if A/R(A) is non- commutative, computing the determinant over A is ModpP-hard. If 
A/R{A) is commutative, then the determinant can be computed in polynomial time. 

Proof of Theorem 5. 5. The algorithm first computes the Wedderburn-Malcev decomposition A = 
B ©' R{A) of the algebra A: a result of de Graaf et al. (Theorem A. 5) shows that such a decom- 
position may be computed efficiently. By Lemma 5.3, we can compute the nilpotency index d of 
the algebra in deterministic polynomial time. We assume that d < n; otherwise, the bruteforce 
algorithm for the determinant has running time N'-^^'^\ 

For any i and j, the (i,j)th entry of the input matrix M can be written uniquely as ruij = 
bij + rij where bij G B and rjj G R] the elements bi^j and rjj are also efficiently computable. 
Now, note that the determinant of the input matrix M can be written as 

det(M) = ^ Sgn(a)(6i^„(i)+ri^^(l))-(62,a(2)+''2,<7(2))---(^n,a(ri)+'^n,a(n)) = X] ^g'^(^) X] ^(^''5') 
o-e5n creSn sc.[n\ 

where t(cr, S) is the product, in increasing order of i, of r,j' for i £ S and bi,j for i ^ S. Note that 
t{a, S) £ R{Ay^^ (we use here the fact that R{A) is an ideal in A) and hence, t{a, S") = if |5| > d. 
Thus, we may only consider S of size strictly less than d. 

We divide the terms t{a,S) based on the r^j that actually appear in t{a,S). Specifically, for 
each 1-1 function / : 5 — )• [n], let t{a, S, f) denote t{a, S) if a\s = f and otherwise. We can write 
the determinant det(M) as 

det(M)= ^ ^ sgn(a)t(a,5,/)= J] det(M(5,/)) 

SC[n]: f:S^[n]: aGSn SC[n]: f:S^[n]: 

\S\<d f 1-1 \S\<d f 1-1 

where the entries m{S, f)ij of M{S, f) are defined as follows: for i £ S, m{S, f)ij = if f{i) ^ j 
and rij otherwise; for i ^ S, m{S, f)ij = bij. We show that for each S and / as above, det(M(S', /)) 
can be computed in time N^^'^\ which will prove the theorem, since there are only vP^'^) of them to 
compute. For the remainder of the proof, we fix some subset 5 C [n] of size t < d and f : S ^ [n] 
that is 1-1. 

Note that the matrix M{S, f) is "almost" a matrix over the commutative subalgebra B of A: 
it contains exactly d entries outside B, one in each row indexed by an element of S. We reduce 
the computation of det(M(5', /)) to the computation of the determinant of a similar matrix over 
a commutative algebra closely related to B. Indeed, let denote B ® B ® ■ ■ ■ ® B (t-|-l 

times). This is a commutative algebra of dimension at most D"^. Furthermore, we see that l®(*+i) 
is the identity element of this algebra. For i G [t] U {0}, we denote by Bi the following subalgebra of 

can easily be seen that each Bi is isomorphic to B by the isomorphism 
4>i: B B^ where (/)i(6) = 1®' O 6 ® 
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Given any z G [n], we denote by Pre(i) the set {i' €z S \ i' < i}. We now construct a new matrix 
M'{S,f) with entries from as follows: 

r if i G 5 and f(i) / j, 

m'{S, f)ij = I ifieS and f{i) = j, 

[ ^e{m{S,f)ij) if i S and ^ = |Pre(i)|. 

In words, to construct M'{S,f), we have replaced each entry in M{S,f) that is in R{A) by the 
identity and each entry bij G -B by the corresponding element in Bi where £ = |Pre(i)|. 

Since M'{S, f) is a matrix with entries from the commutative algebra its determinant 

can be computed in time N'^^^\ Say S = {ii, . . . ,it} and f{ik) = jk for k G [t]. Let {ei, . . . , em} 
be a basis for B. Then, we have 

det(M(5', /)) = sgn(cj) I JJ ^^^^(i) J • u^j, ■ I J]^ ^^^^(i) J • r^^ • • • r^, • f J| 6^ J 

cre5n: \i<n / \ii<i<j2 / \i>it / 



Each product of the form JlieT ^i,o-(i) that appears in the summation above is an element of the 
commutative algebra B and hence can be expanded in the basis of B as follows: 

(m \ / m \ / m \ 

k=l / \k=l / \fc=l / 

^\s=f 

= Y [ Y • • • J ekon^nek, ■ ■ ■ ri.j.Ck, 

ko, — ,kt&[m] \(T£S„:a\s=f / 
~ y ^ '^k ' ^kofil,h^k\ ' ' ' fit,jt^kt (5-1) 

fce[m,]'+i 

where k denotes the tuple {ko,...,kt) and denotes a^^^^ • • • a^*''^. Let us expand 

det(M'(5, /)) similarly. We use e^^^ to denote (/)i{ek)- We have 

det(M'(5,/))= J] sgn(a)(n'^o(6,,.(,))yi^(*+i)-( J] </'i( W)) J " l""^*""'^ " " " l""^*^'^ • ( 11 '^*( W)) 



m \ / m \ / m 

= Y ^g-(-) ( Y 4>f • 1^^*^^^ • Y 4>1^^ • 1^^'-^^^ • • • 1^^*^^^ • Y 4M*^ 

(TeS„: \fc=l / \fe=l / \fc=l 

^\s=f 

= Y "^k- "^k^f^k^ ■ ■ ■ 4t = Y % • efco ^ efc, ® • • • «) ek, 

Thus, we can simply read off the coefficients from det(M'(5', /)), and using Equation (5.1), 
we can compute det(M(S, /)). Since det(M'(5', /)) can be computed in time N'^^'^\ we obtain a 
jV'^('^)-time algorithm to compute det(M(5', /)) and hence for det(M) as well. □ 
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6 Discussion 



Our results show that the basic Godsil-Gutman approach to approximating the permanent, as 
generahzed by Chien et al. [CRS03] runs into many obstacles, since the estimators are not efficiently 
computable. In the case of the quaternions, the result of Chien et al. shows that a suitable 
modification of the basic estimator still gives a relatively good approximation to the permanent. Is 
there such a modification for matrix algebras? 

Our dichotomy theorem in Section 5 used crucially the fact that we worked over a finite field. 
Over the rationals, for example, even the structure of semisimple algebras is fairly complicated, and 
we don't have an exact characterization of when the determinant over such an algebra is efficiently 
computable. Extending our dichotomy theorem to these algebras is an interesting problem. 

Theorem 5.5 shows that even when given the algebra A as input, the determinant remains 
efficiently computable as long as A/R{A) is commutative and A has bounded nilpotence index. How 
close is this to being a characterization of algebras over which the determinant is polynomial-time 
computable (under reasonable complexity assumptions) when the algebra is part of the input? More 
generally, can one come up with suitable conditions on the radical R{A) under which computing 
the determinant over A is hard even when A/R(A) is commutative? 
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A Computing the structure of algebras 

An algebra is simple if it is isomorphic to a matrix algebra (possibly of dimension 1) over a division 
ring containing F. Note that if F is finite, Wedderburn's Little Theorem [Lam91] implies that the 
division ring is a field extension of F and hence, a simple algebra is simply a matrix algebra over a 
field extension of F. An algebra is said to be semisimple if it can be written as the direct sum of 
simple algebras. 

Friedl and Ronyai [FR85] first considered the question of efficiently computing the structural 
elements of an algebra given as input in the form of a multiplication table. That is, the algebra A 
is presented to the algorithm in the form of a basis {ai, . . . , ud} along with a table that lists the 
pairwise products ajCj for i,j € [D]. They proved the following result: 

Theorem A.l ([FR85], Theorem 5.7). There is a deterministic algorithm which, when given a 
finite dimensional algebra A over ¥, computes a basis for the radical R{A). The algorithm runs in 
time poly(A^, log |F|), where N denotes the size of the input. 
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Thus, using the above algorithm, we can obtain in deterministic polynomial time a description 
of A/R{A), the "semisimple" part of A. Friedl and Ronyai [FR85] and Ronyai [R6n87] showed 
respectively that semisimple algebras can further be decomposed into simple algebras, and when 
F is finite, one can find explicit isomorphisms from simple algebras to matrix algebras. We state 
these two results below. 

Theorem A. 2 ([FR85, Theorem 7.8], [R6n87, Theorem 6.2]). There is a deterministic algorithm 
which, when given a finite dimensional semisimple algebra A over a finite field ¥, computes a 
decomposition of A = ^ ■ Ai into simple matrix algebras, and explicit isomorphisms from each Ai 
to a matrix algebra over a field extension of¥. The algorithms run in time poly(A^), where N is 
the size of the input. 

We say that an algebra A is a semidirect sum of subalgebras Bi and B2 \i A = Bi ® B2 as a 
vector space. We will denote this as ^ = i?i ©' -62- The Wedderburn-Malcev theorem tells us that 
any algebra over a finite field or the rationals is a semidirect sum of its radical with a subalgebra. 

Theorem A. 3 ([CR62, Chapter X]). Given any finite dimensional algebra A over a finite field or 
the rationals, there exists a subalgebra B of A such that A = B (B' R{A) . 

A decomposition of the algebra A as given above is called a Wedderburn-Malcev decomposition 
of A. 

Remark A. 4. Note that by definition, given any Wedderburn-Malcev decomposition of A = B (B' 
R{A), the subalgebra B is isomorphic to A/R{A), the "semisimple" part of A. 

The result of de Graaf et al. [GIKR97] shows that given an algebra A and the quotient A/R{A), 
it is possible to obtain a Wedderburn-Malcev decomposition of A = B (B R{A) in deterministic 
polynomial time. We state the result below. 

Theorem A. 5 ([GIKR97, Theorem 3.1]). There is a deterministic polynomial-time algorithm 
which, when given a finite dimensional algebra A over ¥, computes a Wedderburn-Malcev decom- 
position A = B (B' RiA) of the algebra. 

B Computing the determinant over a given algebra 

In this section, we assume that the algebra A is presented to the algorithm in the form of a basis 
{ai, . . . , 0£)} along with a table that lists the pairwise products OiOj for i,j G [D]. We would like to 
efficiently compute the determinant of an input n x n matrix with entries from the algebra A. We 
assume that D < poly(n). (It is easy to see that when D is very large, say n!, even the bruteforce 
algorithm is efficient. So we assume that D is small.) 

Note that under no constraints on the algebra A, this problem is at least as hard as computing 
the determinant over M2(F) and hence the hardness results of Theorem 3.5 apply. We would like 
general conditions on the algebra under which the problem becomes tractable. When F is finite, 
such conditions should ensure that the semisimple part A/R{A) is commutative (or else M2(F) is 
a subalgebra of A). By the Wedderburn-Malcev theorem, this implies that A = B (B' R{A), where 
S is a commutative subalgebra of A. We would like general conditions on R{A) under which the 
determinant is efficiently computable. 

At this point, let us look at an important example that motivates this work. Consider the 
algebra Ud{¥) of dxd upper-triangular matrix algebras over F. Let Dd{¥) denote the (commutative) 
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subalgebra of d x d diagonal matrices and let U'^{¥) be the subalgebra of strictly upper triangular 
matrices (i.e., elements of Ud(¥) that contain only zeroes along the diagonal). It is well-known (see 
[EGH+09, Section 2.5], for example) that R{Ud{¥)) = U'^{¥) and that Ud{¥) = Dd{¥) 0' U'^{¥) is 
a Wedderburn-Malcev decomposition of Ud{¥). Thus, Ud{¥) / R{Ud(¥)) is commutative and hence, 
this is the kind of algebra we would like to consider. It can be shown (and we will see below) 
that when d is constant, the determinant over Ud{¥) can be computed in polynomial time. On 
the other hand, when d is n^^^\ then it is known that the problem is hard; this is implicit in an 
earlier result of Arvind and Srinivasan [ASIO], and can also be proved from Theorem 3.5 above. 
Therefore, any criterion that characterizes algebras w.r.t. the tractability of the determinant over 
them must explain this difference. 

We suspect that this criterion is the nilpotency index of the radical of the algebra in question. 
It is easy to see that the nilpotency index of of U^{¥) is d. In this case, we can use Theorem 5.5, 
which implies the following: 

Corollary B.l. For any constant d and field ¥ that is either finite or the rationals, there is a 
deterministic polynomial-time algorithm running in time poly(A^'^), which when given as input the 
description of an algebra A over ¥ with nilpotence index hounded by d and an nxn matrix M with 
entries from A, computes the determinant of M . 

Thus, the nilpotency index successfully "explains" the tractability of computing the determinant 
over Ud{¥), when d = 0(1). But does the nilpotency index explain every such instance? We suspect 
that the issue of when exactly the determinant becomes intractable is closely related to the cases 
when d becomes large, but we are unable to prove this. 
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